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A formula for the radio ray refraction angle, r, is derived by integration of the approx- 
imate differential equation for the case where the refractivity, (n— 1)X10 6 , decreases expo- 
nentially with height above the surface of a smooth, spherical earth. The solution is in 
terms of the widely tabulated exponential and error functions, and is accurate to within 
4 percent over the useful range of the variables employed. 



1. Background 

In problems where the direction of radio propaga- 
tion is important , e.g., tracking radar, it is convenient 
to be able to compute the refraction of radio rays as 
they pass through the atmosphere. Assuming that 
the refractive index, //, is a function only of height 
above the surface of a smooth, spherical earth, the 
usual integral for the ray refraction angle, usually 
called the "bending", r. is [1, 2] z 



~-r 



cos d {) \(l N/dh]dhXU)- 



sin % ) +~-2(N-N) cos^oXlO" 6 V 

(1) 



where A T is the "target" height to which the bending 
is to be calculated, A T is the radio refractivity, 
(n—l)X10 8 , with N s the surface value, and O is the 
initial elevation angle of the radio ray. The geom- 
etry involved in the problem is shown on figure 1. 
The integral (1) is to be evaluated under the 
assumption that 



N=N s exp {-ch} 



(2) 



where c is a constant determining the gradient of 
N with respect to height; this is the exponential 
model atmosphere recently recommended for inter- 
national tise by the CCIR [3]. 

2. Development 

The integral (1) with N given by eq (2) is put 
in readily integrable form by means of the 
approximation 

iV,XlO- 6 [l-exp {-ch}]~yh (3) 

made where N»—N appears in the denominator of 
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(1). The parameter y is defined by 

N,Xl O-<[l-ex p{-cH}] 
y= - w - 



0<H<h T 



where // is an effective integration height which was 
determined from published tabulations of refraction 
variables in model atmospheres of exponential form 
[4] for 0o = O, Jl,,( l ' s given by 

77=4.75 [exp {— 0.01158ciV,}] fl-exp {-ch T }\. (4) 

With (3), the integral for r in an exponential 
atmosphere is approximated by 



r^iV.XlO^cos 



JO J ai™* 



exp {—ch}dh 



/sin -0,,+2/J T — t cos 2 o V 



(5) 
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Figure 1. Geometry of radio ray refraction. 
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The change of variable: 



where 



2 _ckr sin % 
to- 2 ' 

1— Yr o cos 2 o 



exp 



\-z 2 ]dz. 



results in: 

t ^ N s X 1 - 6 cos 0o V'2 cAr exp { z\ } I 

This can be given in terms of the error function, [5] 

2 



erf (%)=— 

■\Jir Jo 



e-* 2 d£, 



as: 



r^A^XlO^costfoV^exp 



^oj 



[erf{ 2 (A r )}-erf{3„}], (6) 

where -\ 7rr /2, r in kilometers, has been absorbed as a 
factor of 100 in the constant 1CT 4 . 

Equation (6) has the following form for the case of 
total bending, h T ^70 km: 

T(h T ^™)^N s X10- 4 cos6 \rtexp{zl}[l-eYi{z }},(7) 
which reduces to an interestingly short form for d =0, 

3. Accuracy of the Solution 

Equation (6) approximates the true value of r with 
increasing accuracy as the initial elevation angle 
approaches ir/2; this can be shown by substituting 
the asymptotic expansion of the error function [5] in 
eq (7), with the result that 

T(h T ^~)^N s XlO-° COt do [l-£i + £~ a + ■ ■ "I- 

L ^0 4 ^0 OC J 

If z is large enough so that \ z\ may be neglected 



(approximately for O >O.3 radians) then the familiar 
approximation to the total bending at large elevation 
angles results [1, 2], 

r(h T ~^oo)^N s XlO- 6 cot d , O >O.3, 

which is known to be convergent to the true value 
of t as #o approaches t/2. The extension to the case 
of h T <70 km yields similar results. 

The accuracy of eq (6), for 6 <^0.3 radians, has 
been checked against published results obtained for 
the CRPL Exponential Reference Atmospheres [4]. 
The largest errors found, using (4) to obtain H, were 
about 4 percent, and these were for a profile with a 
very strong iV-gradient; the largest errors for profiles 
of average JV-gradient were about 1.5 percent. 

For ease of calculation H may be set constant for 
all profiles and all h T ; a value of H^ 1 km results in 
errors no larger than about 10 percent. 

On the other hand if higher accuracy is desirable 
H may be made a function of initial elevation angle 
also; adding /(^o) to the value of H as given by (4), 
where 

/(0o) ^ i-L 2 4fl ' °° in radians ' 

will limit the errors to about 1 percent for the CRPL 
exponential profiles. 

In summary, eq (6) represents a relatively accurate 
and concise formulation for the ray bending in an 
exponential atmosphere. It has the advantage of 
being in terms of elementary, well-tabulated func- 
tions, and so can be easily programed for solution on 
digital or analogue computers, and in addition the 
solution is asymptotic to the correct values for large 
initial elevation angles. 
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